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Summary 


Hereditary xerocytosis is thought to be a rare genetic condition characterized by red blood cell 
(RBC) dehydration with mild hemolysis. RBC dehydration is linked to reduced Plasmodium 
infection in vitro; however, the role of RBC dehydration in protection against malaria in vivo is 
unknown. Most cases of hereditary xerocytosis are associated with gain-of-function mutations in 
PIEZOI, a mechanically activated ion channel. We engineered a mouse model of hereditary 
xerocytosis and show that Plasmodium infection fails to cause experimental cerebral malaria in 
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these mice due to the action of Piezol in RBCs and in T cells. Remarkably, we identified a novel 
human gain-of-function P/EZO/ allele, E756del, present in a third of African population. RBCs 
from individuals carrying this allele are dehydrated and display reduced Plasmodium infection in 
vitro. The existence of a gain-of-function PIEZO1 at such high frequencies is surprising, and 
suggests an association with malaria resistance. 


Graphical Abstract 


A gain of function mutation in the mechanically activated channel PIEZO1 is associated with 
resistance to the malaria parasite Plasmodium falciparum 


— $ 
C C S 
Wild type Piezo1 GOF Wild type E756del carrier 
(Piezo1 GOF) 
Normal Ded Normal x me 
RBC BC RBC BC 


LE E & 


*Plasmodium *Plasmodium *Plasmodium *Plasmodium 


IL LOW 
P. falciparum P. falciparum 
infection rate infection rate 


Cerebral malaria protected eC— 2 


Introduction 


PIEZOs are non-selective cation channels that sense mechanical stimuli in many 
multicellular organisms (Coste et al., 2010; Ranade et al., 2015). PIEZO1 is essential for 
mechanotransduction in vascular development, blood pressure regulation, and red blood cell 
(RBC) volume control, among other roles (Li et al., 2014; Ranade et al., 2014; Retailleau et 
al., 2015; Wang et al., 2016; Cahalan et al., 2015). The related PIEZO2 is the principal 
mechanosensor for touch and proprioception (Ranade et al., 2014; Woo et al., 2014; Woo et 
al., 2015; Chesler et al., 2016). Human genetic studies have highlighted the significance of 
PIEZO1 in human development and physiology. Patients with loss-of-function mutations in 
PIEZO1 suffer from persistent lymphedema caused by congenital lymphatic dysplasia 
(Lukacs et al., 2015). PIEZO/ mutations are also linked to hereditary xerocytosis, also 
known as dehydrated hereditary stomatocytosis, (Zarychanski et al., 2012; Albuisson et al., 
2013; Bae et al., 2013). Hereditary xerocytosis is a dominantly inherited blood disorder 
characterized by RBC dehydration causing reduced RBC osmotic fragility, and is associated 
with mild or asymptomatic hemolysis (Delaunay, 2004). This disorder is considered to be 
rare, and found mostly in the Caucasian population (Archer et al., 2014; Glogowska et al., 
2017). Complications include splenomegaly, resulting from increased RBC trapping in the 
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Results 


spleen, as well as iron overload due to unknown mechanisms (Archer et al., 2014). 19 
different point mutations in PIEZOI have been described to cause hereditary xerocytosis 
(Murthy et al., 2017). Some of these mutations have been electrophysiologically analyzed, 
and show slower inactivation kinetics compared to wildtype PIEZO1 channels. The slower 
inactivation translates to more ions passing through PIEZO ion channels, and thus these 
mutations are considered gain-of-function. Consistently, Piezol deficiency in RBCs in mice 
causes overhydration (Cahalan et al., 2015). Beyond RBCs, studies in mice have suggested 
wide-ranging functions of PIEZO1 in various biological processes; whether hereditary 
xerocytosis is associated with other conditions beyond RBC pathology is currently not fully 
understood. 


Plasmodium, the causative parasite for malaria, has exerted strong selective pressures on the 
human genome (Kwiatkowski, 2005). This is demonstrated by severe genetic conditions, 
such as sickle cell disease, that persist in human populations from malaria-endemic areas 
because the underlying genetic variants confer resistance to Plasmodium infection (Hedrick, 
2004; Feng et al., 2004). The scope of RBC disorders that might contribute to Plasmodium 
resistance, however, has not been fully explored. Interestingly, dehydrated RBCs (including 
those from hereditary xerocytosis patients) show delayed infection rates to Plasmodium in 
Vitro, suggesting a potential protective mechanism against infections from this parasite 
(Tiffert et al., 2005). The effects of dehydrated RBC on Plasmodium infection in vivo, 
however, remain unknown. Since overactive PIEZO1 causes dehydrated RBCs in hereditary 
xerocytosis patients, we reasoned that mice carrying a gain-of-function Piezo/ allele could 
offer a suitable model to investigate the effect of Plasmodium infection in vivo (de Oca et 
al., 2013). 


Piezo? gain-of-function mice recapitulate human hereditary xerocytosis phenotypes 


To test whether gain-of-function Piezo/ expression causes xerocytosis-like phenotypes in 
mice, and to elucidate the role of xerocytosis in Plasmodium infection in vivo, we 
engineered mice that conditionally express a human-equivalent hereditary xerocytosis 
mutation (Figure 1). Specifically, R2456H is a xerocytosis mutation in human P/EZO/ that 
displays significantly longer channel inactivation time (t) (Albuisson et al., 2013). We 
verified that the equivalent mouse Piezo] point mutation (R2482H), when overexpressed in 
HEK cells that lack endogenous PIEZO/ (PIEZOIKO HEK) (Dubin et al., 2017), showed 
slower channel inactivation (Figure 1A). Since residue 2482 resides in the last coding exon 
(51), we designed the knock-in construct by flanking exons 45-51 with /oxP sites, followed 
by a copy of the region containing exons 45-51 with a mutation that would replace R with H 
at residue 2482 (Figure 1B). We named this conditional allele Piezo/™. In cells that express 
Cre recombinase, the wild type exon will be replaced by the modified exon, allowing tissue- 
specific control of gain-of-function Piezo] expression. 


We generated a constitutive gain-of-function Piezo/ mouse line by crossing mice 
homozygous for the mutant allele (Prezo/€***) with cmv-cre mice that expressed a Cre 
driver ubiquitously (Schwenk et al., 1995). We also generated a hematopoietic lineage- 
specific gain-of-function Piezo/ mouse line (PiezoJGOF blood ) using vav1-cre (de Boer et 
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al., 2003). To evaluate the expression of the gain-of-function allele, we sequenced the last 
exon of Piezo! cDNA from whole blood of homozygous PiezoIGOF »!°°4, and observed the 
expected nucleotide change c.GG7742.7743AC (Figure S1A). In addition, we found that 
Piezol transcript levels in whole blood from both homozygous and heterozygous 

Piezo 1GOF Plo*d mice were similar to levels observed in wild type mice, demonstrating that 
the genetic manipulation did not alter Piezo/ expression levels (Figure S1B). We also found 
that both constitutive (P/ezo/GOF constitutive) and blood cell-specific (Piezo/GOF blood) 
transgenic mice (heterozygous and homozygous) were born at the expected Mendelian ratio 
and appeared to develop normally. 


We found that RBCs from both homozygous and heterozygous Piezo /GOF blood mice 
showed reduced osmotic fragility, as shown by a left-shifted curve in a hypotonicity- 
dependent hemolysis challenge (Figure 1C-D). This demonstrates that RBCs from gain-of- 
function Piezo/ mice are more resistant to lysis in response to hypotonic solutions compared 
to wild type, a defining feature for hereditary xerocytosis (Archer et al., 2014). Piezo/GOF 
blood mice also displayed hematological properties similar to mild anemia, indicated by a 
lower hemoglobin level and increased reticulocyte number, as is the case for individuals with 
hereditary xerocytosis (Table S1). Those patients also have increased Mean Corpuscular 
Volume, which is a measure of RBC volume, and increased Mean Corpuscular Hemoglobin, 
which indicates average hemoglobin mass per RBC (Zarychanski et al., 2012; Albuisson et 
al., 2013; Bae et al., 2013; Archer et al., 2014). We found that shifts in these two values in 
gain-of-function Piezo] mice were similar to those observed in hereditary xerocytosis 
patients. Mean Cell Hemoglobin Concentration, in contrast — which is expected to be 
elevated in dehydrated RBCs - was not significantly increased in homozygous Piezo/GOF 
blood mice. Importantly, however, we found that RBCs from these mice were nevertheless 
dehydrated, as they showed reduced osmotic fragility (Figure 1C—D). In addition, we used 
scanning electron microscopy and found the presence of RBCs with deformed and 
dehydrated shapes from heterozygous Piezo /GOFl?9d mice, which is another clinical 
feature often observed in patients (Figure 1E). One of the predominant features of hereditary 
xerocytosis is splenomegaly. We found that both homozygous and heterozygous Piezo/GOF 
blood mice had significantly larger spleens (1.04+0.03cm2, 0.74+0.02cm2, respectively) 
compared to wild types (0.4240.02cm2, n= 4 animals per genotype, Student's t-test, 
compared to wild type, p«1x107^) (Figure 1F and Figure S1C). Together, our data show that 
gain-of-function Piezo] mice display hallmark clinical features observed in human 
hereditary xerocytosis patients, including RBC dehydration, mild anemia, and splenomegaly. 


Gain-of-function Piezo1 mice have reduced growth rate of Plasmodium blood stages and 
protect against experimental cerebral malaria 


To evaluate the connection between Piezo/, RBC dehydration, and protection against 
malaria, we infected gain-of-function Piezo/ mice with a GFP-expressing reference line of 
the ANKA strain of rodent malaria parasite Plasmodium berghei (Franke-Fayard et al., 
2004). We chose P berghei ANKA since this parasite is a well-established model to analyze 
the course of infections 7n vivo, and to investigate experimental cerebral malaria in mice 
(Franke-Fayard et al., 2004; de Souza et al., 2010; Hunt et al., 2010). We found that wild 
type mice died between day 6 and 8, consistent with previous findings (Franke-Fayard et al., 
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2004; de Oca et al., 2013) (Figure 2A). In contrast, we observed that the homozygous and 
heterozygous Piezo /GOF constitutive mice survived as long as 24 and 19 days, respectively 
(Figure 2A). Importantly, the post-infection survival rates of Piezo JGOF Plo9d mice were 
indistinguishable from Piezo/GOF °"Sttutive mice, indicating that induction of gain-of- 
function Piezo/ in hematopoietic lineages was sufficient to extend post-infection survival 
(Figure 2A). 


Next, we analyzed the course of infection in wild type and gain-of-function Prezo/ mice to 
test if the expression of the mutant Piezo/ allele affects Plasmodium growth rate in RBCs, as 
suggested by previous in vitro experiments (Tiffert et al., 2005). We measured the 
percentage of RBCs that were GFP-positive (parasitemia) by flow cytometry. During the 
first week of infection (Phase 1, Figure 2B), we found that parasitemia reached 6-12% in 
wild type mice at the time of death; however, both Piezo/GOF constitutive and Piezo IGOF 
blood mice had significantly lower parasitemia (on day 6, 5.14+0.42% for constitutive mice 
and 5.20+0.34% for blood cell-specific mice, p<0.05 compared to wild type, 8.53+1.65%, 
Student's t-test). These findings suggest that expression of a gain-of-function Prezo/ allele 
in blood cells reduce parasite growth rate of blood stages (Figure 2B). Unlike wild type 
animals, which all died at the end of Phase 1, gain-of-function Prezo/ mice then entered a 
second phase of infection (Phase 2; day 7 to day 18, Figure 2C). We found that during this 
phase, they exhibited a steady increase in parasitemia, eventually leading to severe 
hyperparasitemia of up to 70% of infected RBCs (Figure 2B and 2C). These data suggest 
that gain-of-function Piezo/ expression can dramatically modify the course of Plasmodium 
infection in vivo, leading to enhanced survival, despite high end-stage levels of parasitemia 
(Figure 2C). 


A prominent feature of experimental cerebral malaria in the P berghei ANKA/C57BI/6 
infection model is the breakdown of blood-brain barrier (Nacer et al., 2014). We injected 
Evans blue dye into mice and studied blood-brain barrier compromise. As expected, we 
observed blue dye leakage into brain parenchyma in all wild type mice at day 6 after 
infection (n=8) (Figure 2D, left), indicating blood-brain barrier breakdown. Remarkably, we 
did not detect Evans blue leakage in the brains of Piezo/GOF blood (n=7) even at day 18 
when they were about to die (Figure 2D, right). To quantitatively evaluate blood-brain 
barrier disruption, we measured the optical density of Evans blue dyes extracted from 
infected brains (n = 5 per genotype) (Ferreira et al., 2011). We observed a significant 
reduction in brain Evans blue contents in infected Piezo/GOF^l?9d compared to wild type 
mice (Figure 2E). In addition, we evaluated experimental cerebral malaria by measuring 
brain water content that reflects the severity of brain edema caused by cerebral 
complications (Hunt et al., 2014). Wildtype mice had increased brain water content after 
infection compared to Piezo/GOF Plo9d mice (Figure 2F). Thus, our data show that gain-of- 
function Piezo/-expressing mice are protected against experimental cerebral malaria. 
However, these mice eventually died, probably due to severe anemia, as they showed 
reduced HGB levels (2.85+2.5 g/dL, n=3, in PiezoIGOF »!°4 mice 18 days after infection) 
compared to uninfected Piezo /GOF Plo?d mice (14.02 0.16 g/dL, n=5, p«0.002) (Phillips 
and Pasvol, 1992). Together, our results suggest that gain-of-function Piezo/ expression 
reduces Plasmodium growth rate of blood-stage infection in vivo and can protect mice from 
the development of cerebral complications. The reduced Plasmodium infection rate of 
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dehydrated RBCs observed 7n vitro (Tiffert et al., 2005) can explain the reduced parasite 
growth rate of blood stage observed in gain-of-function Piezo/ mice during phase 1; 
however, a connection between dehydrated RBCs and protection from experimental cerebral 
malaria was novel and unexpected. 


Red blood cell dehydration is responsible for reduced parasite growth and partially 
responsible for protection against cerebral malaria in gain-of-function Piezo1 mice 


To address whether decreased parasite growth rate and prevention of experimental cerebral 
malaria in the gain-of-function Piezo/ mice were due to RBC dehydration, we genetically 
rescued RBC dehydration in gain-of-function Piezo/ mice, and assessed P berghei infection. 
We took advantage of the fact that PIEZO1-induced RBC dehydration requires the activity 
of KCa3.1, a calcium dependent potassium channel (also known as Gardos channel). 
Activation of KCa3.1 drives potassium and water out of RBCs in response to increased 
intracellular calcium, thereby causing dehydration (Maher and Kuchel, 2003; Cahalan et al., 
2015). We crossed the gain-of-function Piezo/ mice to KCa3./ knockout mice. As expected, 
Piezo 1GOF blood; KCa3.1~~ mice had osmotic fragility similar to wild type mice, 
demonstrating that RBC dehydration was corrected by removing KCa3.1 channel activity 
(Figure 3A and Figure S2). After P berghei ANKA infection, PiezoIGOF 5lood/ KCa3. 17^ 
mice survived significantly longer than wild type (but shorter than Piezo/GOF Plood), This 
resulted in an intermediate survival curve of Piezo/GOF blood/ KCa3. 1 mice (p<0.0001, 
compared to wild type and Piezo 1GOF Plo?) (Figure 3B). This suggests that correction of 
RBC dehydration fails to reverse survival rate to a level that is similar to wild type, 
indicating that RBC dehydration is not completely responsible for the increased survival rate 
in the gain-of-function Piezo/ mice. 


We also found that Piezo/GOF ^lood/ C53, |- mice had a parasite growth rate of blood- 
stage that was indistinguishable from that of wild type during the first week of infection, 
suggesting that RBC dehydration was responsible for the reduced parasite growth observed 
during Phase 1 in gain-of-function Piezo/ mice (Figure 3C). Importantly, KCa3. / knockout 
mice in wild type Piezo/ background did not show changes in parasitemia, suggesting that 
the absence of KCa3./ per se did not influence RBC infection (Figure 3C, gray). Finally, 
quantitative measurements of both Evans blue and brain water content in infected brains 
showed that PiezoJGOF ^lood/ KCa3. [- mice experienced an intermediate level of cerebral 
complications between wild type and Piezo/GOF Plo9d mice (Figure 3D-3F). Together, our 
data from PiezoIGOF blood; KCa3.1~“ genetic experiments suggest that 1) RBC 
dehydration is completely responsible for the reduced parasite growth rate (Phase 1); and 2) 
RBC dehydration is a major contributing factor for the absence of experimental cerebral 
malaria (Phase 2), but that other mechanisms may be involved. 


Gain-of-function Piezo1 expression in red blood cells and T cells contributes to protection 
against cerebral malaria 


The incomplete protection from experimental cerebral malaria in Píezo/GOF blood; 
KCa3. [—— mice (despite normal parasite growth rate) suggests the existence of other 
mechanisms that affect cerebral complication in gain-of-function Piezo/ mice. Previous 
work has shown that processes critical for the development of cerebral malaria in both 
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humans and rodents involve both RBCs and immune cells (Baptista et al., 2010; Nacer et al., 
2014; Dunst et al., 2017). To directly address the cell autonomous function of gain-of- 
function Piezo/ allele in these cells, we induced expression of gain-of-function Piezo/ 
mutation in different blood cell types and tested survival rate, parasite growth rate, and 
experimental cerebral malaria. 


First, we generated RBC-specific gain-of-function Piezo/ mice (Piezo IGOF RBO) with 
EpoR-cre (Heinrich et al., 2004). We verified the efficiency and specificity of EpoR-cre 
expression by measuring RBC osmotic fragility for Piezo/GOF RBC mice. We found that 
these mice had reduced RBC fragility, similar to Píezo/GOF Plo9d mice, suggesting that 
EpoR-cre was efficiently inducing recombinase activity in most RBCs (Figure 4A and 
Figure S3A). Also, gain-of-function Piezo/ mRNA was not present in immune cells (CD4+ 
and CD8+ T cells) from Piezo/GOF RBC mice. This is an important control, as we address 
the role of gain-of-function Piezo/ expression in T cells separately (see below) (Figure 
S4B). We found that infection of PiezoIGOF RBC mice with P berghei caused a survival rate 
indistinguishable from Piezo/GOF *l?9d mice (Figure 4B). Furthermore, Piezo/GOF RBC 
mice had a parasitemia curve indistinguishable from Piezo/GOF Plo9d mice, and did not 
develop experimental cerebral malaria (Figure 4C and Figure D). These data suggest that the 
expression of gain-of-function Piezo/ in RBCs is sufficient to cause reduced Plasmodium 
growth rates and to protect mice from the development of cerebral complications. 


Parasite-specific CD8+ cells are essential in causing Plasmodium-induced cerebral 
complications (Yañez et al., 1996; Belnoue et al., 2002; Howland et al., 2015). We induced 
gain-of-function Piezo/ expression in peripheral CD4+ and CD8- T cells by using hCD2- 
cre (Vacchio et al., 2014). We tested the specificity of hCD2-cre by measuring RBC osmotic 
fragility and showed that Prezo/°/*; hCD2-cre (PiezoIGOF Tells) mice had normal RBC 
fragility, confirming that hCD2-cre did not induce gain-of-function Piezo/ expression in 
RBCs (Figure 4A and Figure S3A). Furthermore, we evaluated the efficiency of hCD2-cre in 
targeting CD4+ and CD8- T cells by sequencing the cDNA made by those cells from 
homozygous PiezoJGOF T-e!l5 mice, and found that gain-of-function Piezo? mRNA was the 
only Piezo/ transcript expressed in the targeted cells (Figure S3B). 


We found that after P berghei infection, PiezoIGOF Tels mice survived significantly 
longer than wild type mice (p«0.01), but not as long as Piezo/GOF Plood or Piezo /GOF RBC 
mice (p«0.01), suggesting that expression of gain-of-function Prezo/ in T-cells provided 
partial protection (Figure 4B). Furthermore, we found that parasitemia in Piezo/GOF Tells 
mice was identical to that of wild type mice during the first 7 days after infection, before it 
continued climbing until the end of the infection (Figure 4C, compare blue and black). This 
suggested that gain-of-function Prezo/ expression in CD4/8- T cells did not alter parasite 
growth rate of blood stage compared to wild type mice (Phase 1). Intriguingly, despite wild 
type-like parasite growth rates during the first 7 days, Piezo/GOF TeS mice displayed 
attenuated experimental cerebral malaria during Phase 2 (Figure 4D and 4E). Also, 
PiezoIGOF T-*!l5 mice had an intermediate level of cerebral complications between wild 
type and PiezoIGOF RBC mice (Figure 4F). These data demonstrate that gain-of-function 
Piezol expression in T cells can provide partial survival advantage by attenuating the 
disruption of the blood-brain barrier seen in experimental cerebral malaria. 
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Finally, we tested whether gain-of-function Piezo/ expression in macrophages can affect 
Plasmodium infection, since macrophages have been shown to be important for both 
protection and pathology in malaria (Chua et al., 2013). We expressed gain-of-function 
Piezol specifically in macrophages using LysM-cre (Clausen et al., 1999)). PiezoIGOF 
macrophage mice displayed survival rate (Figure 4B) and parasitemia curves similar to wild 
type littermates (Figure 4C). These results indicate that macrophages are unlikely to play an 
essential role in reducing parasite growth rate and protection against experimental cerebral 
malaria in xerocytosis mice. Together, our data suggest that RBCs play a major role in gain- 
of-function Piezo/-mediated protection against Plasmodium infection and cerebral malaria; 
however, T cells also appear to be involved in protection against cerebral complications. 


Identification of a common human PIEZO7 gain-of-function mutation in African 
populations under positive selection 


The role of gain-of-function Prezo/ in rodent malaria described here raises a conundrum: If 
Piezol mutations are protective against Plasmodium infection, why then is hereditary 
xerocytosis not commonly observed in individuals from Africa, where malaria is highly 
prevalent? To investigate whether common P/EZO/ gain-of-function mutations are present 
in African populations, we took a comparative genomics approach to look for possible 
PIEZO1/ gain-of-function alleles, and catalogued nonsynonymous (missense) single 
nucleotide polymorphisms (SNPs) and in-frame insertions/deletions (indels) in PIEZO/. To 
maximize the likelihood of finding gain-of-function mutations, we: 1) performed our search 
using the Exome Aggregation Consortium data (ExAC) (Lek et al., Exome Aggregation 
Consortium, 2016); 2) picked PIEZO/ SNPs or indels with allele frequencies above 0.5%; 
and 3) picked mutations that were more than 5-fold enriched in African populations, as 
compared to people of non-African descent. Using these criteria, we found 21 mutations 
consisting of 19 SNPs and 2 indels (Table S2). 


To test for potential functional effects of the various mutations, we next performed a large- 
scale calcium-imaging assay by screening the 21 candidate mutations for increased response 
to various concentrations of the PIEZOI agonist Yodal (Syeda et al., 2015). We found that 
two of the 21 mutations lead to increased PIEZO1 responses in this screen: amino acid 
substitution A1988V (SNP) and indel E756del (3 nucleotide deletion) (Figure 5A-B). We 
found that the A1988V mutation only has an allelic frequency of 0.8% in the African 
population (inset in Figure 5A). In contrast, the E756del mutation has an allelic frequency of 
18% in individuals of African descent (3% in Europeans), and therefore present in at least 1 
copy in about a third of African population (inset in Figure 5B). 


To test if these mutations lead to gain-of-function PIEZOI channel kinetics, we recorded 
mechanically activated currents and found that PIEZO1 variants containing A1988V or 
E756del mutations was activated normally by mechanical force, but had significantly longer 
inactivation time constants (t) compared to wild type (p«0.0001). This is similar to R2456H, 
a gain-of-function allele that has the longest inactivation time among all hereditary 
xerocytosis mutations (Albuisson et al., 2013) (Figure 5C-D), and the equivalent of this 
allele was used to create our gain-of-function Piezo/ mice (Figure 1). These data show that 
gain-of-function PJEZO! mutations with similar ion channel activities to those causing 
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hereditary xerocytosis in Caucasian families (Albuisson et al., 2013) can be found in 
individuals of African descent. At least one of these, E756del, is present in one third of 
African individuals, suggesting a potential connection between PIEZO1, hereditary 
xerocytosis, and malaria. 


We focused on the more abundant allele, E756del. We hypothesized that this allele may be 
under positive selection in African populations, where malaria is endemic. To test this 
hypothesis, we assessed three main signatures of selection, commonly found in allelic 
variants under positive selection (Sabeti et al., 2006): 1) population differentiation of the 
allele observed between African and non-African populations, as measured by 7st; 2) 
whether a variant is in linkage-disequilibrium with nearby SNPs creating a long-range 
haplotype block, which is commonly observed in more recent («25,000 years) selective 
sweeps; and 3) whether the allele is derived (7.e. non-ancestral), since such new alleles 
typically have low population frequencies, unless under selection. 


We looked at the frequency of the E756del allele in the populations present in the 1000 
Genomes catalogue (1000 Genomes Project Consortium, 2015), and found that it is present 
at high allelic frequency (9% to 23%) in all populations of African descent, including 
African Americans (allelic frequency 14%), but not in individuals of non-African ancestry 
(allelic frequency «146, Figure 6A). The observed genotype frequencies at this locus are in 
Hardy-Weinberg equilibrium f =0.201, p = 0.654), and therefore segregating as expected 
in a randomly mating population. Next, we investigated population differentiation across the 
1000 Genome populations. We calculated Foz values at all PIEZO/ missense mutation loci 
between individuals of African and non-African descent, and found that the populations 
were most differentiated using the E756del allele (Fs for E756del = 0.32, Fer of all other 
PIEZO! missense alleles = 0—0.26, Figure 6B). This finding is consistent with E756del 
being under positive selection in populations where malaria is endemic. 


We next investigated the regions surrounding the E756del locus, but did not observe any 
SNPs in significant linkage disequilibrium with E756del. The lack of an observed long 
haplotype flanking this allele makes it harder to conclusively provide proof of positive 
selection of the E756del variant (Vitti et al., 2013). The lack of linkage disequilibrium, 
however, could also be because this allele might have been subject to selection on standing 
variation (że. not as the result of a selective sweep (Sabeti et al., 2006)), or because the 
selective pressure on this locus is relatively old (» 25,000 years). Even though Plasmodium is 
an ancient parasite (Loy et al., 2017), the former is still a likely explanation because the 
expansion of P falciparum and subsequent impacts on human selection likely began in the 
last 10,000 years (Joy et al., 2003). 


To assess whether the E756del variant is derived (i.e. is a new allele that occurred in 
Africans) or ancestral, we investigated the architecture of the PJEZO/ locus in the archaic 
humans and non-human primates. There is low amino acid sequence homology near the 
E756del locus between humans and non-human primates (Figure S4A), thus we could not 
investigate pre-human ancestry. We found, however, that both Neanderthals and Denisovans 
had the wild type E756 in their PIEZO/ genes (Figure 6C, and Figure S4B). This finding 
shows that the PIEZO/ E756del gain-of-function allele is derived in individuals of African 
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descent, again consistent with being under positive selection (Sabeti et al., 2006). Combined, 
our analyses show that the PJEZO/ gain-of-function mutation E756del is a high-frequency 
(present in 1/3" of African population) derived allele that is highly differentiated in 
populations where malaria is endemic. These findings are highly suggestive of the E756del 
genetic variant being under positive selection in populations of African descent (Sabeti et al., 
2006), presumably because of its likely role as a malaria-protective allele. 


Red blood cells from E756del African American carriers are dehydrated and cause reduced 
infection by Plasmodium falciparum in vitro 


We acquired blood samples from healthy volunteer African American donors and tested 
whether E756del causes xerocytosis-like RBC dehydration, and importantly, whether it 
confers attenuation of infection against P falciparum in vitro. We obtained 25 whole blood 
samples and used white blood cells to sequence the exon containing E756del. We found that 
nine (3696) African American donors were heterozygous for E756del (none were 
homozygous) (Figure 55A). We also screened all 25 donors for other known common 
mutations that affect RBC morphology and could potentially influence susceptibility to 
Plasmodium infection. Our sequencing results showed that all 25 donors were free of HbS, 
HbC and HbE mutations in the B-globin chains that cause hemoglobinopathy (Figure S5B). 
In addition, we showed that none of the 25 donors had the variations that cause a.- 
thalassemia (Figure S5C), another condition associated with RBC abnormality and 
Plasmodium infection (Chong et al., 2000). 


Next, we imaged RBCs with scanning electron microscopy from three carriers and showed 
that all had RBCs with echinocyte and stomatocyte morphologies, which is a characteristic 
of hereditary xerocytosis RBCs (Figure 7A). Remarkably, we also found that RBCs from all 
9 donors with the E756del mutation were dehydrated as assayed by osmotic fragility test 
(Figure 7B-C), similar to RBCs from known xerocytosis patients (Delaunay, 2004; Archer 
et al., 2014). Next, we infected both control and E756del carrier RBCs with P falciparum in 
Vitro. Parasitemia was significantly lower for E756del carriers relative to non-carriers, 
measured by both Giemsa and SYBR green staining methods (Figure 7D, E) (Johnson et al., 
2007). Together, our data demonstrates that E756del is a common PIEZO/ gain-of-function 
mutation in African populations, causing RBC dehydration in otherwise healthy African 
Americans, and is likely under positive selection, due to its ability to confer reduced 
susceptibility of RBCs to P falciparum infection. 


Discussion 


Gain-of-function Piezo1 expression in blood cells provides protection against 
Plasmodium-induced cerebral complications in vivo 


Dehydrated RBCS, including those from hereditary xerocytosis patients, show slower 
infection rates to P falciparum in vitro (Tiffert et al., 2005). However, this mechanism of 
protection has never been tested in vivo. To address these issues, we engineered a gain-of- 
function Piezo/ mouse that recapitulated most features of hereditary xerocytosis. 
Remarkably, gain-of-function Piezo/ mutation induced in different types of blood cells 
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caused dramatic shifts in survival rates in response to P berghei infection, caused by reduced 
parasite growth rate of blood stage as well as protection from experimental cerebral malaria. 


Our mouse genetic data suggests that gain-of-function Piezo/-induced RBC dehydration is a 
major determinant in the protection against cerebral complications of malaria. Several other 
genetic mutations that affect RBC morphology are associated with resistance to malaria in 
human populations (Hedrick, 2004; Feng et al., 2004), and some of these mutations also 
cause RBC dehydration, such as sickle cell disease (Brugnara, 1995). Similar experiments 
can be performed in the future to evaluate the potential contribution of RBC dehydration to 
malaria resistance in the genetic disorders mentioned above. Another important next step is 
to determine the molecular mechanisms responsible for RBC dehydration-dependent 
attenuation of Plasmodium infection. 


In addition to RBC dehydration, we discovered an unexpected function of gain-of-function 
PIEZO1 in immune cells during Plasmodium infection. T cells play both pathogenic and 
protective roles in human malaria, as well as in murine malaria models (Hafalla et al., 2006; 
Ewer et al., 2013). T cells experience diverse mechanical stimuli during development and 
function, but the role of mechanosensitive ion channels in immune cells are poorly 
understood (Huse, 2017). It is possible that overactive PIEZO1 channels alter T cell 
developmental programs and/or modulate their activity when encountering parasites. It will 
be of interest to use both gain-of-function and loss-of-function Piezol mice to explore the 
role of this ion channel in T cells. 


The discovery of gain-of-function PIEZO!1 allele present in one third of African population 


The discovery of gain-of-function PIEZO/ E756del in African populations with a high allele 
frequency of ~18% (such that an estimated third of African people carry this mutation as 
heterozygotes) is quite surprising. Our findings dramatically redefine the epidemiology of 
this disorder: hereditary xerocytosis-like condition is much more common than previously 
anticipated. Thus, E756del provides a unique opportunity to evaluate the association 
between gain-of-function PIEZO1, RBC dehydration, and malaria in endemic regions. 


Despite the experimental evidence above, PIEZO/ locus was not detected as a strong 
candidate by recent genome-wide association studies (GWAS) that aimed to identify genetic 
loci for severe malaria resistance (Malaria Genomic Epidemiology Network and Malaria 
Genomic Epidemiology Network, 2014). This is potentially due to GWAS limitations and 
the complexity of this particular genetic locus. GWAS samples have high levels of genetic 
diversity and are underrepresented in reference panels of genetic variation (Leffler et al., 
Malaria Genomic Epidemiology Network, 2017). Also, GWAS studies mainly use SNPs to 
determine association and this would be challenging to evaluate more complex loci without 
genetic imputation method. E756del is in such a locus with multiple short tandem repeats 
(Figure 6C), so that imputation of this mutation into current GWAS datasets is not 
straightforward. In this regard, our experimental data provide promising clues for association 
analysis: sequencing this particular locus in endemic population can determine whether 
E756del is associated with protection against severe malaria. 
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E756del provides an opportunity to evaluate the role of over-active mechanotransduction 
in human health 


Does E756del allele cause hereditary xerocytosis and other disorders? We readily identified 
E756del carriers from self-reported healthy African American blood donors. Whether 
E756del carriers have anemia or splenomegaly is not known to date. A full clinical 
evaluation of individuals carrying this allele will be of high interest to assess how overactive 
PIEZO1 influences xerocytosis-related phenotypes, as well as other conditions. For example, 
analysis of loss-of-function Piezo/ mice has demonstrated a critical role of this ion channel 
in cardiovascular function (Retailleau et al., 2015; Wang et al., 2016; Rode et al., 2017). 
Therefore, it will be of interest to assess the role of overactive PIEZOI channel in 
hypertension, which has high incidence in African Americans (Kaplan, 1994). We expect 
that a complete clinical characterization of individuals with the E756del allele will shed 
further light on the range of phenotypes that are associated with PIEZO/, including anemia, 
splenomegaly, autoimmune diseases, various aspects of cardiovascular function, as well as in 
indications not previously associated with PIEZOI. 


STAR Methods 


Contact for Reagents and Resource Sharing 


Further information and requests for resources and reagents should be directed to and will be 
fulfilled by the Lead Contact, Ardem Patapoutian (ardem @scripps.edu). 


Experimental Model and Subject Details 


All animal procedures were approved by the Institutional Animal Care and Use Committees 
of The Scripps Research Institute (TSRI). 


Mice—PiezoIGOF Id, Piezo IGOF constitutive. Piezo IGOF Teels, Piezo IGOF macrophage, 
and PiezoIGOF RBC mice were generated by breeding Piezo/cx/cx with vav1-cre (The 
Jackson Laboratory, stock# 018968) and cmv-cre (The Jackson Laboratory stock# 006054), 
hCD2-cre (The Jackson Laboratory stock# 027406), LysM-cre (The Jackson Laboratory 
stock# 004781), and EpoR-cre (a gift from Dr. Klingmuller group at Max-Planck-Institute 
für Immunbiologie, Freiburg, Germany). KCa3. 1 —/— mice were ordered from The Jackson 
Laboratory (stock# 018826). Gain-of-function Piezo/ mice were generated and maintained 
on C57BL/6 background. All animals were backcrossed at least 10 generations to C57BL/6. 
The mice were housed in a 12hr light/dark cycle (light from 6am to 6pm) in a temperature- 
controlled room (24 degree) with free access to food and water. The ages and sexes of mice 
are indicated in the following method section. Littermates were used for experiments. 


Cell lines and cell culture—P/EZO/KO HEK cells were grown in Dulbecco's modified 
Eagle's medium containing 4.5 mg/ml glucose, 1096 fetal bovine serum, 1x antibiotics/ 
antimycotics. 


Human blood samples—The collection of human whole blood was approved by TSRI 
normal blood donor service. Fresh whole blood from TSRI normal blood donor service was 
drawn and kept at ambient temperature (in heparin-coated containers), followed by osmotic 
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fragility test and P. falciparum infection experiments (see below) on the same day. Whole 
blood from Biological Specialty Corporation was delivered by air at ambient temperature the 
next day after collection for further experiments. 


Method details 


P. berghei infections and parasitemia measurement by flow cytometry—Donor 
mice (C57BL/6) were intraperitoneally injected with 50—200ul of erythrocytes parasitized 
with either P berghei (ANKA) GFPcon 259cl2 (provided by California Institute for 
Biomedical Research, La Jolla, USA) or P berghei (ANKA) mCherry-hsp70-Luc-eefla (line 
1868 from Leiden Malaria Research Group, the Netherlands). We used P berghei (ANKA) 
mCherry-hsp70-Luc-eefla when infecting blood cell-specific GOF Piezo] mice because one 
of the Cre driver (EpoR-cre) has EGFP expression in RBCs so that P berghei (ANKA) 
GFPcon 259cl2 cannot be used. Blood was collected by cardiac puncture from infected 
donors when the parasitemia reached 4—6% (see below). Parasitized erythrocytes were 
washed with sterile saline three times at 1000xg for 3min and diluted to 5x10? infected 
cells/ml as working solution. 200ul working solution was intravenously injected into the 
experimental mice for analysis. For GFP-fluorescence based parasitemia measurement, 1.5ul 
tail blood was collected from infected mice in 180pl Dulbecco's Phosphate Buffer Solution 
with 2% fetal bovine serum, on 96-well plates. The cytometry was performed on NovoCyte 
Flow Cytometer system (ACEA Biosciences, San Diego, CA) following manufacturer’ s 
instructions. Briefly, erythrocytes were selected on size for analysis by gating on forward/ 
side-light scatter. Excitation of erythrocytes was performed with a laser at a wavelength of 
488 nm and emission of the green fluorescence, or a wavelength of 587nm and emission of 
the mCherry was detected using a filter of 530 nm. By gating the uninfected erythrocytes 
and the GFP-positive infected erythrocytes parasitemia was calculated as the percentage of 
infected cells. 


Blood-brain barrier and experimental cerebral malaria assay—When P berghei 
infected displayed either experimental cerebral malaria (including ataxia, convulsion, 
paralysis and/or coma) or severe anemia (immobility and pale blood color), 2% Evans Blue 
(Sigma-Aldrich, dissolved in sterile saline) was intravenously injected into P berghei 
infected at 5ml/kg body weight. After 45-60 min, euthanized mice were transcardially 
perfused with Phosphate Buffer Solution and 4% paraformaldehyde before brains were 
dissected. To quantify the Evan blue (EB) contents, the infected brains (n = 5 animals per 
genotype) were cut into small pieces and incubated in 1ml formamide at room temperature 
for 36 hours, followed by measuring the optical density at 620nm by Cytation3 (BioTek, 
Winooski, VT). The concentration of EB dyes was calculated from a standard curve with the 
equation (Y = 0.03263*X + 0.03413, where Y= reading, X= EB concentration). To quantify 
the brain water content, the infected mice (n = 5 animals per genotype) were sacrificed and 
brains were dissected. Brain tissues were weighted (wet weight), then dehydrated at 56?C. 
The sample was re-weighted 48hr later to obtain a dry weight. The percentage of water was 
calculated by: BWC = [(wet weight — dry weight) / wet weight] x 100%. Note there are 
cases that single brains were cut into half, with each half was used for both Evans blue and 
water content calculations. 
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Scanning Electron Microscope—Samples of mouse blood were added to ice cold 
buffered saline (10mM NaCl, 155mM KCl, 10mM glucose and 1 mM magnesium chloride) 
before being fixed in 2.5% glutaraldehyde in 0.1 M cacodylate buffer on ice. Aliquots of the 
fixed cells were placed on 12mm coverslips previously coated in polylysine for 30 mins. 
Following a buffer wash and postfixation in buffered 1% osmium tetroxide, the samples 
were extensively washed in distilled water. The samples were dehydrated in graded ethanol 
series followed by processing in a critical point drier (Tousimis autosamdri 815). The 
coverslips were then mounted onto SEM stubs with carbon tape and sputter coated with 
Iridium (EMS model 150T S) for subsequent examination and documentation on a Hitachi 
S-4800 SEM (Hitachi High Technologies America Inc., Pleasanton CA) at 5kV. 


For human RBCs, samples were fixed at room temperature with 2.5% glutaraldehyde in 
0.1M cacodylate buffer followed by 1% osmium tetroxide, and washed in water. A 10 ul 
drop of suspension was loaded on the sample carrier and imaged in a FEI Quanta200 FEG 
microscope in ESEM mode using the gaseous secondary electron detector. The stage was 
set-up at 2?C, the acceleration voltage was 15kV and the working distance 10 mm. Water 
was then progressively removed by cycles of decreasing pressure/injection of water, until 
reaching equilibrium at the dew point. The minimal final pressure in the chamber was 350 
Pa. Pictures were taken with a dwell time of 6 psec. 


Osmotic fragility test and hematology test—Blood (for both mouse and human) was 
diluted at 1:8 into isotonic saline (0.9% NaCl) containing 2 mM HEPES, pH 7.4. 10 ul of the 
diluted blood was pipetted into each well (in a row) on a 96-well round-bottom plate. 
Separate rows were used for separate samples. 225 ul tonicity solutions made from saline 
solutions at concentration of 0, 20, 25, 30, 35, 40, 45, 50, 55, 60, 80, and 100%. The plate 
was incubated for 5min at room temperature followed by 5 min centrifuge at 1000xg. 150 ul 
supernatant was transferred to 96-well flat bottom plate for absorbance reading at 540nm 
using Cytation3 (BioTek, Winooski, VT). The data was analyzed using 4-parameter 
sigmoidal nonlinear regression. Hematological properties from mice were analyzed by 
hematology analyzer at Ruggeri lab, Dept. of Molecular and Experimental Medicine, The 
Scripps Research Institute, La Jolla, USA) 


Gain-of-function Piezo1 mice generation—The targeting strategy was based on the 
NCBI transcript NM. 001037298.1. Wildtype exons 45-51, including the complete 3' 
untranslated region (UTR) were flanked by loxP sites. An additional polyadenylation signal 
(nucleotide sequence of the synthetic polyA: 
gagctccctggcgaattcggtaccaataaaagagctttattttcatgatctgtetettgetttttetetecgecgcg) was inserted 
between the 3^ UTR and the distal loxP site in order to prevent downstream transcription of 
the mutated exon 51 in the conditional allele. The size of the loxP-flanked region is 2.8 kb. 
The exons 45-51, including the splice acceptor site of intron 44 were duplicated and inserted 
downstream of the distal loxP site. The R2482H mutation was introduced into the duplicated 
exon 51. Positive selection markers were flanked by FRT (Neomycin resistance - NeoR) and 
F3 (Puromycin resistance- PuroR) sites and inserted into intron 44 and downstream of the 
synthetic polyA, respectively. The targeting vector was generated using BAC clones from 
the C57BL/6J RPCIB-731 BAC library and transfected into the TaconicArtemis C57BL/6N 
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Tac ES cell line (Taconic, Hudson, NY). Homologous recombinant clones were isolated 
using double positive (NeoR and PuroR) selection. The conditional knock-in allele was 
obtained after Flp-mediated removal of the selection markers. The constitutive knock-in 
allele was obtained after Cre-mediated deletion of wildtype exons 45-51 and the synthetic 
polyA sequences 


Mechanical stimulation—For whole-cell recordings, mechanical stimulation was 
achieved using a fire-polished glass pipette (tip diameter 3—4 um) positioned at an angle of 
80? to the recorded ce Ils. Downward movement of the probe towards the cell was driven by 
a Clampex controlled piezo-electric crystal microstage (E625 LVPZT Controller/Amplifier; 
PhysikInstrumente). The probe had a velocity of 1 um.ms~! during the ramp segment of the 
command for forward motion and the stimulus was applied for 150 ms. To assess the 
mechanical responses of a cell, the probe was first placed as close to the cell as possible (this 
distance could vary from cell to cell). Then, a series of mechanical steps in 1 um increments 
was applied every 10s, which allowed full recovery of mechanically activated (MA) 
currents. Threshold was calculated by subtracting the distance at which the probe first 
touched the cell surface from the minimal distance at which mechanically activated currents 
were evoked. Mechanically activated inward currents were recorded at a holding potential of 
—80 mV. The inactivation kinetics at a holding potential of -80 mV of traces of currents 
reaching at least 75 96 of the maximal amplitude of current elicited per cell were fitted with 
mono-exponential equation (or in some case bi-exponential equation for the rapidly-adapting 
currents, accordingly to previous reports (Albuisson et al., 2013) and using the fast time 
constant giving a value of inactivation time (x) per responsive cell used for analysis. Channel 
kinetic properties between WT and mutant PIEZO1 were compared using Student's t-test. 


Cell culture and transient transfection—P/EZO/KO HEK cells were grown in 
Dulbecco’s modified Eagle's medium containing 4.5 mg/ml glucose, 10% fetal bovine 
serum, 1x antibiotics/antimycotics. Cells were plated in 6-well plates and transfected using 
lipofectamines 2000 (Invitrogen by ThermoFisher Scientific, Carlsbad, CA), according to 
the manufacturer's instruction. Human P/EZO/ or mouse Piezol mutations fused to IRES- 
td Tomato was transfected at 1.4ug per well (6-well plate) for fluorescent imaging plate 
reader (see below). To measure calcium signals, ultrasensitive sensor GCaMP6 (Chen et al., 
2013) were transfected at 0.6ug per well (6-well plate). Cells were incubated for 2 days 
before electrophysiology experiments or fluorescent imaging plate reader. 


Fluorescent imaging plate reader (384-well format)—A fter transfection, the cells 
were dissociated from 6-well plates and re-seeded into 384-well plate, at 12,000 cells per 
well. The plate was incubated for 2 days then washed with assay buffer (1x HBSS, 10 mM 
HEPES, pH7.4) using a ELx405 CW plate washer (BioTek, Winooski, VT). Fluorescence 
was monitored on a fluorescent imaging plate reader (FLIPR) Tetra. A 10-mM stock 
solution of Yodal in dimethyl sulfoxide (DMSO) was used resulting in a maximum of 196 
DMSO in the assay. 10 uM Yodal was used in initial screens for searching gain-of-function 
mutations (compared to wild type). Positive hits were then validated by using a series of 
Yodal concentrations. Concentration-response curves were fitted using a sigmoidal dose— 
response at variable slope (GraphPad Prism, La Jolla, CA). 
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Real time quantitative PCR— Total RNA was isolated from mouse whole blood by 
Quick-RNA Whole Blood (Zymo Research, Irvine, CA). 500 ng total RNA was used to 
generate Ist strand cDNA using the Quantitect reverse transcription kit (Qiagen). Real time 
PCR assays were set up using GoTaq qPCR Master Mix (Promega, Madison WI). The 
reaction was run in the ABI 7900HT fast real time system using 1 ul of the cDNA in a 20 ul 
reaction according to the manufacturer's instructions in triplicates. Primers were designed 
for target gene (mPiezo/) and reference gene (Gapdh). See Key Resource Table for primer 
information. Calibrations and normalizations were done using the 2-AACT method, where 
AACT = ((CT(target gene) -CT (reference gene)) — (CT (calibrator) — CT (reference gene)). 


CD4+ and CD8+ T cell isolation—MoJoSort Mouse CD4 and CD8 Nanobeads kits 
(BioLegend, San Diego, CA) were used for magnet-based cell separation of CD4/8+ T cells. 
We performed the procedures based on instructions provided by the kit manual. 


Population genetic analysis—We obtained minor allele and genotype frequencies from 
the Exome Aggregation Consortium (ExAC) and the 1000 Genomes Project. 2504 genomes 
were analyzed, 661 from African and 1843 from non-African ancestries. Wright's fixation 
index (Fey), a measure of population differentiation, was calculated as follows: 


2 
o 
s 


F on = — = 
ST p-p) 


where pis the mean allele frequency and o? is the allele frequency variance between the 


populations. The 1000 Genomes browser (http://phase3browser. l000genomes.org) was used 
to determine that no alleles were in linkage disequilibrium with E756del (estimated r? values 
were <0.05). Hardy-Weinberg equilibrium was estimated using the classical binomial 
expansion to determine the expected genotype frequencies and x? tests. 


PIEZO! sequences were obtained from the following sources: modern humans (Homo 
sapiens, Genbank NG_042229.1), neandertals (Homo neanderthalensis, Neandertal Ensembl 
ENSG00000103335), and Denisovans (Homo sapiens ssp. Denisova, previously generated 
reads (Meyer et al., 2012; Reich et al., 2010) aligning to humans [GRCh37/hg19] using the 
UCSC Genome Browser, (Kent et al., 2002)). All non-human primate amino acid sequences 
were obtained from GenBank. The sequences were aligned in Geneious using MAFFT 
(Katoh, 2002). 


Genotyping in African American blood donors 


Genotyping E756del carriers by sequencing: 25 whole blood samples (5—10ml) were 


collected from Normal Blood Donor Service (The Scripps Research Institute, La Jolla, CA) 
and Biological Specialty Corporation (Colmar, PA), approved by institutional regulations. 
200ul of whole blood samples were used for genomic DNA isolation by QIAamp DNA 
Blood Mini Kit (Qiagen, Germany). A -200bp PCR amplicon that contained E756 locus 
was generated for sequencing E756del or wild type allele (forward primer: 
5’CAGGCAGGATGCAGTGAGTG3 ’, reverse primer: 
5’GGACATGGCACAGCAGACTG3’. Reverse primer was used for sequencing). 
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Screening for hemoglobin mutations by sequencing: PCR amplicons that contained 


potential sickle cell mutation sites were generated (Figure S5). Forward primers: 
5' AGAAGAGCCAAGGACAGGTA3'; reverse primers: 
5’TTGCAGCCTCACCTTCTTTC3’. Reverse primer was used for sequencing. 


Screening for a-thalassemia by multiplex PCR: Each 50 uL reaction contained 20 
mmol/L Tris-HCl pH 8.4, 50 mmol/L KCI, 1.5 mmol/L MgCl2, 1 mol/L betaine (SIGMA, 
St. Louis, MO), 0.2 uL of each primer (see below) 0.2 mmol/L of each dNTP, 2.5 units of 
polymerase, and 50-100 ng of genomic DNA. Reactions were carried out on a thermal 
cycler, with an initial 5-minute denaturation at 95°C, 30 cycles of 97°C for 45 seconds, 60°C 
for 1 minute 15 seconds, 72?C for 2 minutes 30 seconds, a nd a final extension at 72?C for 5 
minutes Primers: 


a2/3.7-F CCCCTCGCCAAGTCCACCC, 3.7/20.5-R AAAGCACTCTAGGGTCCAGCG, 
a2-R AGACCAGGAAGGGCCGGTG, 4.2-F GGTTTACCCATGTGGTGCCTC, 4.2-R 
CCCGTTGGATCTTCTCATTTCCC, SEA-F CGATCTGGGCTCTGTGTTCTC, SEA-R 
AGCCCACGTTGTGTTCATGGC, FIL-F TGCAAATATGTTTCTCTCATTCTGTG, FIL-R 
ATAACCTTTATCTGCCACATGTAGC, 20.5-F GCCCAACATCCGGAGTACATG, MED-F 
TACCCTTTGCAAGCACACGTAC, and MED-R TCAATCTCCGACAGCTCCGAC 


P. falciparum culture—P falciparum Dd2 strain parasites were cultured under standard 
conditions (Trager and Jensen, 1976), using RPMI media supplemented with 0.05 mg/ml 
gentamycin, 0.014 mg/ml hypoxanthine (prepared fresh), 38.4 mM HEPES, 0.2% Sodium 
Bicarbonate, 3.4 mM Sodium Hydroxide, 0.05% O+ Human Serum (Denatured at 56°C for 
40 min; Interstate Blood Bank, Memphis, TN) and 0.0025% Albumax). Human O+ whole 
blood was obtained from TSRI Normal blood donor service (La Jolla, CA). Leukocyte-free 
erythrocytes are stored at 50% hematocrit in RPMI-1640 screening media (as above, but 
without O+ human serum and with 2x albumax concentration) at 4°C for one to three weeks 
before experimental use. Cultures were monitored every one to two days via Giemsa-stained 
thin smears. 


Parasitemia Determination—Asynchronous P falciparum parasites (Dd2 strain) were 
cultured in standard conditions (as described above), then synchronized twice via sorbitol 
(Lambros and Vanderberg, 1979) and grown to 7% parasitemia at the late trophozoite/ 
schizont stage. Patient blood was obtained from TSRI Normal blood donor service, washed 
and centrifuged three times (at 800 x g for 5 min at 4°C) with RPMI only and once with 
complete media (as above) (Elias and Greene, 1979), with any visible buffy coat being 
removed after each spin. All blood samples were given a numerical designation and allele 
status was not determined until after all data collection was completed. All genotypes were 
blinded from experimenters. (Basic Malaria Microscopy, 2010) 


Parasite growth was determined in two independent ways: absolute parasitemia 
determination via thin-blood smear after Giemsa staining, and inferred parasite growth via 
the SYBR Green I-based fluorescence assay (Johnson et al., 2007). For the absolute 
parasitemia determination by thin blood smear, parasites were established, starting from the 
7% parasitemia cultures above and diluted using the corresponding patient RBCs, in a 10 
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mL culture at 5% hematocrit and 0.5% parasitemia (the same as for the Giemsa absolute 
parasitemia determination assay). Cultures were then grown for 4 days and parasitemia 
measurements were taken every 2 days. For the SYBR green I inferred parasitemia 
determination assay, when parasite burden is estimated based upon DNA content, parasites 
were cultured in 100 uL volumes in 96-well plates at 5% hematocrit and 0.5% parasitemia 
with at least 5 replicates per time point. Parasites were plated on three 96-well black plates 
with clear bottoms (Fisher Scientific). (One surrogate plate was used for measurement every 
two days, with DNA content determined by SYBR Green I incorporation of lysed parasites). 
Relative parasitemia was determined by fluorescence measurement, background was 
determined using uninfected RBCs and subtracted, then relative parasite burden was 
determined via normalization against a known O+ WT blood sample. In both cases all 
measurements were taken for all samples, genotypes were then assigned to numbered patient 
samples, wild type vs heterozygote samples were averaged at each time point, and average 
parasitemia values were compared. 


Quantification and statistical analysis 


Statistical analysis—AlII of the data are presented as the mean + SEM or SD and 
represent at least 3 independent experiments. Statistical analysis, significance level and n 
values are described in the Figure or Figure legends. For mouse experiments n = number of 
animals and at least n = 4 were used. For human blood experiment, n = number of individual 
blood samples. For comparison, we performed two-tailed Student’s t-test, where p < 0.05 is 
considered statistically significant. For all datasets, we used Prism7 to perform the statistical 


analysis. 


Supplementary Material 


Refer to Web version on PubMed Central for supplementary material. 
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Highlights 
Expression of a gain-of-function Prezo/ allele models hereditary xerocytosis 
in mice 
Mice expressing gain-of-function Piezo/ allele are protected from cerebral 


malaria 


A third of African population carry a PIEZO/ gain-of-function allele 
(E756del) 


RBCs from E756del carriers are dehydrated and show reduced susceptibility 
to Plasmodium 
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Mouse model for human xerocytosis. (A) Representative traces of mechanically-activated 
(MA) inward currents for wild type and Piezol R2482H. ****p « 0.001. (B) Strategy for 
generating knock-in mouse. (C) Osmotic fragility test for RBCs. (D) Quantification for 


osmotic fragility. Relative tonicity at which 50% RBCs are lysed (half hemolysis) was 


calculated for each curve. ****p « 0.001. (E) Scanning electron microcopy images. 


Heterozygous Piezo/GOF blood RBCs showed signs of stomatocytes. (F) Splenomegaly in 


gain-of-function Piezo] mice (Fig S1C). Scale bar: 10m. Data are presented as means + 
SEM. See also Figure S1 and Table S1 
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Plasmodium infection in gain-of-function Piezo/ mice. (A) Survival curves for gain-of- 
function Piezo/ mice after P berghei infected RBCs (B) and (C) Parasitemia recorded by 
flow cytometry for phase 1 (first 7 days) and phase 1 and 2 together (24) days respectively. 
(D) Intact blood-brain barrier in infected gain-of-function Prezo/ mice. (E) Quantification of 
blood-brain barrier disruption. (F) Brain water content in infected brains. (*p«0.05, 
***p«0.01, and ***p<0.001). Scale bar: 5mm. Data are presented as means + SEM. 
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Role of RBC dehydration in Plasmodium infection in mice. (A) Deletion of KCa3. / in 
heterozygous PíezoJGOF Plo9d mice (orange) restored RBC dehydration in heterozygous 
Piezo 1GOF blood (green). Piezo/GOF bloody KCa3, 17^ mice had a similar survival curve to 
wild type. (B) Post infection survival rate of Piezo/GOF ^loed/ KCa3. [- mice (orange) is 
intermediate between wild type (black) and heterozygous Píezo/GOF Plo?d mice (green). 
p<0.0001, (Mantel-Cox tests). (C) Both Piezo/GOF Plood/ &C43, 17^ and KCa3. [^- mice 
had same parasitemia as wild type, with significantly higher than heterozygous Prezo/GOF 
blood mice. 

(D) Breakdown of blood-brain barrier in Piezo/GOF Plood/ KCa3.1~“~ mice 13 days after 
infection. (E) Quantification for blood-brain barrier disruption. (F) Brain water content in 
infected brains. (*p«0.05, **p«0.01, and ***p«0.001). Scale bar: 5mm. Data are presented 
as means + SEM. See also Figure S2. 
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Role of gain-of-function Piezo/ expression in RBCs and T cells during Plasmodium 
infection in mice. (A) RBC osmotic fragility for different gain-of-function Piezo/ mice. (B) 
Mice with gain-of-function Piezo/ in different blood cells had distinct survival rates after 
infection. Piezo /GOF RBC (red) had survival rate similar to pan-blood cell-specific mice 
(Piezo IGOF Plo9d (green)), p>0.05. Macrophage-specific gain-of-function mice (PiezoIGOF 
macrophage) had same survival rate as wild type, p>0.05. PiezoJGOF Tells had a survival rate 
greater than wild type (p<0.01) and less than Piezo /GOF Plo9d mice (p«0.01). Mantel-Cox 
tests. (C) Parasitemia recorded by flow-cytometry for gain-of-function mice. *p«0.05, 
**p«0.01, and ***p«0.001, student t-test. (D) Blood-brain barrier compromise in T cell- and 
RBC-specific gain-of-function mice after infection. (E) Quantification for blood-brain 
barrier disruption. (F) Brain water content in infected brains. Scale bar: 5mm. Data are 


presented as means + SEM. See also Figure S3. 
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Identification of gain-of-function P/EZO/ mutations in African populations. (A) and (B) 
Yodal induced Intracellular calcium signals in PIEZO/KO HEK cells overexpressing 
A1988V and E756del cDNA (* p<0.05). Allele frequency for both mutations is shown in the 
insets. (C) Representative traces of mechanically activated (MA) inward currents for wild 
type and mutated cDNA. R2456H, A1988V, and E756del mutations. (D) Quantification for 
inactivation time (t). n = number of cells. (***p < 0.001, **p < 0.01). Data are presented as 


means + SEM 
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PIEZO!1 codon position (minus strand) 
748 750 752 754 756 758 760 762 
Neanderthal CAGCAGGAGGAGGAGGAGGMGGAGGAGRACT (sequence data unavailable) 
Denisovan CAGCAGGAGGAAGAGGAGGAGGAGGAGGACTCCAGGGACGAGGGG 
Modem human CAGCAGGAGGAGGAGGAGGAGGAGGAGGACTCCAGGGACGAGGGG 
Human E756del CAGCAGGAGGAGGAGGAGGAGGAGGA— CTCCAGGGACGAGGGG 
Amino acid Q Q E E|E EE E E DIS RD EG 
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0.0: Deletion derived in 
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Figure 6. 


Population genetics of PIEZO/ gain-of-function E756del allele common in populations of 
African descent. (A) Human population demographics for E756 indels. E756 deletion 
(TCC/-) exists at high frequencies in all populations of African descents (purple). A minor 
allele, E756 insertion (TCC/TCCTCC) was also discovered (coral). (B) Differentiation (FS) 
between populations of African and non-African ancestry at each loci for all PIEZO/ 
missense mutations. Alleles are colored by whether the minor allele frequency (MAF) was 
highest in African (red) or non-African (black) populations, or were similar (grey). (C) A 
nucleotide alignment of modern and pre-modern (Neandertal and Denisovan) human 
PIEZO1 minus strand sequences around the E756del allele showing the codon positions. 
The TCC deletion (GGA on minus strand) spans two codons, but only deletes E756 while 
shifting nucleotides to leave D757 intact. Individual Neandertal and Denisovan reads used to 
create this alignment and comparisons to non-human primate PIEZO1 amino acid sequences 
are shown in Figure S4. 
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Characterization of RBCs from E756del carriers for xerocytosis-like phenotypes, and P 
falciparum infection. (A) SEM images. Three individual E756del carriers have RBCs with 
echinocytes (white arrowhead) and stomatocytes (yellow arrowhead), magnified in lower 


panels. Scale bar for upper panels, 10um; for lower panels, Sum. (B) and (C) Osmotic 
fragility test. RBCs from E756del heterozygous carriers (n=9) had a left-shifted curve (blue) 
compared to controls (n=16) (black), as quantified in (C) **p «0.01. (D) and (E) P 
falciparum infection into RBCs from E756del carriers. Giemsa staining (D) and SYBR 
Green labeling of parasite DNA inside RBCs (E) (**p « 0.01, *p < 0.05). Statistics: 
student's t-test for each time point. Scale bar: 8mm for upper panels in A; 5mm for lower 


panels in A. Data are presented as means + SEM. See also Figure S5. 
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KEY RESOURCES TABLE 


REAGENT or RESOURCE SOURCE IDENTIFIER 


Healthy whole blood samples TSRI normal blood donor service (La Jolla, https://nbds.scripps.edu/ 
CA) 


biological specialty corporation (PA) http://www.biospecialty.com/ 
Chemicals, Peptides, and Recombinant Proteins 
Yodal (50mg stock) (Tocris) Catalog#5586/50 


GoTaq qPCR Master Mix Promega, Madison Catalog#A6002 
e EN NNI 
Experimental Models: Cell Lines 


Human Piezo/KO HEK cells Our own lab Dubin et al., 2017 
ooo Qe 


Experimental Models: Organisms/Strains 


Mouse: B6.129P2- Lyz2m!(crollfoj] The Jackson Lab Stock# 004781 


Mouse: Eport™!(EGFP/cre)Uk Dr. Klingmuller group at Max-Planck- N/A 
Institute für Immunbiologie, Freiburg, 
Germany 


Mouse: Piezolcx/cx Taconic Biosciences Customized 


Plasmodium Berghei: P berghei (ANKA) GFPcon California Institute for Biomedical Research, N/A 
259c12 La Jolla, USA 


Plasmodium Berghei: P berghei (ANKA) mCherry- Leiden Malaria Research Group, the line 1868 
hsp70-Luc-eefla Netherlands 


Plasmodium falciparum Elizabeth Winzeler, university of California, N/A 
San Diego 


Oligonucleotides 
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REAGENT or RESOURCE 


GGTTTACCCATGTGGTGCCTC 
CCCGTTGGATCTTCTCATTTCCC 
CGATCTGGGCTCTGTGTTCTC 
AGCCCACGTTGTGTTCATGGC 
TGCAAATATGTTTCTCTCATTCTGTG 
ATAACCTTTATCTGCCACATGTAGC 
GCCCAACATCCGGAGTACATG 
TACCCTTTGCAAGCACACGTAC 
TCAATCTCCGACAGCTCCGAC 


Recombinant DNA 


Geneious 
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This paper 


GraphPad software http://www.graphpad.com 
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